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Context. High-mass microquasars consist of a massive star and a compact object, the latter producing jets that will 
interact with the stellar wind. The evolution of the jets, and ultimately their radiative outcome, could depend strongly 
on the inhomogeneity of the wind, which calls for a detailed study. 

Aims. The hydrodynamics of the interaction between a jet and a clumpy wind is studied, focusing on the global wind- 
and single clump-jet interplay. 

Methods. We have performed, using the code Ratpenat, three-dimensional numerical simulations of a clumpy wind 
interacting with a mildly relativistic jet, and of individual clumps penetrating into a jet. 

Results. For typical wind and jet velocities, filling factors of about <; 0.1 are already enough for the wind to be considered 
as clumpy. An inhomogeneous wind makes the jet more unstable when crossing the system. Kinetic luminosities ~ 
10"^^ erg/s allow the jet to reach the borders of a compact binary with an O star, as in the smooth wind case, although 
with a substantially higher degree of disruption. When able to enter into the jet, clumps are compressed and heated 
during a time of about their size divided by the sound speed in the shocked clump. Then, clumps quickly disrupt, 
mass-loading and slowing down the jet. 

Conclusions. We conclude that moderate wind dumpiness makes already a strong difference with the homogeneous wind 
case, enhancing jet disruption, mass-loading, bending, and likely energy dissipation in the form of emission. All this 
can have observational consequences at high-energies and also in the large scale radio jets. 

Key words. Hydrodynamics - Radiation mechanisms: general - X-rays: binaries - ISM: jets and outflows - Stars: winds, 
outflows - Shock waves 



1. Introduction 

Microquasars are binary systems hosting a star and an 
accreting black hole or neutron star. Matter from the 
star is transferred to the compact object, part of it 
being launched through magnetocentrifugal forces (e.g., 
Blandford fc Znajek|1977l|Blandford fc Payne|1982l|Barkov 



& Khangulyan 20111. This triggers the formation of bi- 
polar jets, which gen e rate non-the rmal radio emission (e.g. 



Part of the energy carried by the jet can be dissipated 
in the form of magnetic reconnection, recolhmation and in- 
ternal shocks, shear layers in the jet walls, and turbulence. 
Part of the dissipated energy can go to non-thermal parti- 
cles, generating low- and high-energy emission via different 
mechanisms, synchrotron from radio to X-rays, and inverse 
Compton (IC) and hadron-related proc esses up to gamma 
rays (see, e.g. 



Mirabel et al. 1999 Ribo 2005), and are thought to be 



the location from where t he gamma rays observed in some 



erences therein 
As shown W 



Bosch-Ramon & Khangulyan 2009i and ref- 



sources are emitted (e.g., Albert et al. 
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). Jets could 



magnetically dominated 



Perucho fc Bosch-Ramonl ^008^ (PB08 



hereafter) and Perucho et al. i PoiOaj ) (PBKIO hereafter), in 



at 



hydrodynamical processes occurring at higher jet height 
would accelerate the flow, e fliciently converting magn etic 



energy into kinetic one (e.g. Komissarov et al. 2007 1. At 
the scales of the binary system (~ 10°iisch, where i^sch is 
the Schwarzschild radius), the jet is likely to be already a 
hydrodynamical (HD) flow. We focus here on the persis- 
tent jets thought to be present during the low-hard state 
of microquasars, although some considerations for transient 
ejecta, associated to low-h ard to high-soft state tran sitions. 



are done below (see, e.g., Fender et al. 
views on microquasar states) . 
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microquasars hosting an OB star (high-mass microquasars; 
hereafter HMMQ) the jet may be strongly influenced by the 
stellar wind. The one-side impact of the wind on the (pre- 
sumably) already HD jet leads to strong and asymmetric 
recollimation shocks, bending and different types of insta- 
bilities: the recollimation shocks seem suitable candidates 
for particle acceleration and non-thermal emission; bending 
may be noticeable in radio at milliarcsecond scales; instabil- 
ities may destroy the jet flow even within the binary system. 
For typical wind and jet velocities, say Ww ~ 2 x 10^ and 
Wj ~ 10^*^ cm/s, respectively, studies show that for compact 
binaries and jet-to- wind momentum flux ratios ^0.1 the 
jet can be already disrupted (PB08, PBKIO). This number 
is rather constraining, since only a few HMMQ might be 
above this threshold. This could be the reason for the low 
number of HMMQ detected, as suggested in PBKIO. In any 
case, even if not destroyed within the binary system, jets 
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can suffer strong perturbations with both dynamical and 
radiative consequences. 

Previous work in HMMQ wind-jet interactions was done 
under the assumption that the wind is homogene ous, but in 
fact stellar w inds are thought to be clumpy (e.g., Owocki & 
Cohen|2006 ( Moffat 2008 ) . For this reason, it has been pro- 
posed that wind dumpin ess should be taken into account 
when studying HMMQ ( .Owocki et al ||2009[ lAraudo et al. 
20091 -ABR09 hereafter- rlRomero et al.|2010l Araudo et al. 



20111. An important parameter that determines the inho- 
mogeneity of the wind is the wind filling factor /, which de- 
termines the wind volume fraction with higher density. For 
a significant departure from homogeneity, the intraclump 
medium mass and momentum fluxes will be negligible and 
only clumps will have a dynamical impact on the jet. Since 
the interaction between a HD jet and a clumpy wind has not 
been studied in detail, we have carried out 3-dimensional 
(3D) simulations of this scenario. Simulations have been 
done for two different jet powers and jet-to- wind momen- 
tum ratio, Lj = 3 X 10^^^ - lO^^^ erg/s and w 0.03 - 0.08, re- 
spectively, to explore what could be the transition between 
jet destruction and long-term coUimation. Another simu- 
lation has focused on the evolution of individual clumps 
injected in the jet at different heights. Unlike in PBKIO, in 
which the simulation started with the jet being injected at 
its base, here the jet is conical and crossing the whole grid, 
and the clumpy wind is injected from one of the jet sides. 

The first goal of this work is the study of the hydrody- 
namical evolution of a jet when the wind interacting with 
it is clumpy. The second goal is to quantify for which values 
of the clump and jet parameters, dumpiness becomes a rel- 
evant factor. It is also interesting to study the evolution of 
a clump under the impact of a microquasar jet. The results 
can also be used to refine radiation models or interpret ra- 
dio observations, although this will be treated qualitatively. 
The paper is organized as follows: in Section 2 the scenario 
studied here is briefly in trod uced; in Sect. 3 the simula- 
tions are described (Sect. 3.1); results are shown in Sect.llj 
finally, in Sect.jSj the results are discu ssed in the context of 
jet propagation in HM MQ (Sect. 5.1), individual clump-jet 



interactions (Sect. [572|), and their implications for the non- 
thermal emission (? 



5.3 ). Throughout the paper, we will 



use cgs units. 



2. Physical scenario 

The scenario studied here consists of a jet crossing the bi- 
nary system in a HMMQ. The jet starts close to the com- 
pact object, which is located at a distance oid = 2x 10^^ cm 
(following PBKIO) from the massive star, and is perpendic- 
ular to the orbital plane. The jet is initially conical, with 
a radius to height ratio of r; = Rj/z — 0.1. The scenario 
is similar to that studied in PBKIO. However, unlike in 
that work, the stellar wind is assumed here to be inho- 
mogeneous, with a filling factor / ~ 0.1. A sketch of the 
considered scenario is presented in Fig. [ij 

The inhomogeneities or clumps are modelled as gaus- 
sians with ct = i?c = 3 x 10^° cm In reality, the size, 
mass, and velocity of clu mps in stellar winds may follow 
complex distributions (e.g., Moffat|2008[ ). However, we have 
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Fig. 1. Sketch of the scenario simulated in this work. The 
green dashed line confines the simulated region. 



adopted average and representative values for these quan- 
tities for simplicity. Bigger and thus less numerous clumps 
than assumed here would probably have a stronger impact 
on the jet dynamics, unless they were so few that interac- 
tions were rare. The latter may be the case, provided that 
wind mass seems to concentrate in the small clumps . These 
small clumps may also be denser than big ones (Moffat 



2008). On the other hand, if i?c < 3 x 10^° cm, then the 
wind could be effectively considered as homogeneous (as 
in PBKIO). The reason is that clumps cover inside the jet 
a fraction of its radius: x — Rc/R] ^ Vw Rc ("" '^jPc/Sij)^^^ 
(for a Newtonian jet), where Pc = Al/'iirfcPv^ is the clump 
density. This estimate is based on the clump disruption 
time, expected to be slightly longer than the clump shock 
crossing time: ^ Rc/cc, where ^ {2 Lj/nVjPcRfy^^ is 
the shocked clump sound speed. 

Another relevant timescale is the clump acceleration 
time along the jet, i.e., the time needed to accelerate the 
clump material up to a speed ~ Cc, which is ~ td, and 
because of quick expansion, it takes only several times 
this value for the clump material to reach ~ wj (see, 
e.g.,|Blandford fc Koenigl|[T979l [Klein et al.|[T994l ABR09; 



Pittard et al.||2010| [Barkov et a 



2010 



for shocked clump 
For stellar mass-loss rates 
/s, d ^ 2 X 10^^ cm, and 
10^0 cm) (//O.I)- 



-1/2 



^ The effective size of such clumps, if they were spherical and 
homogeneous, would be ~ 4 x 10^° cm. 



evolution in different contexts) 
M - lO-<^Af0/yr, Lj - lO^^ erg/ 
~ 10^° cm/s, X is 0.1 (i?c/3 
i.e., clumps with i?c <C 3 x 10^° (//O.l) cm are destroyed in 
the external jet layers. In this small clump case the wind-jet 
contact discontinuity should probably develop a turbulent 
shear-layer faster than in the homogeneous case. Due to 
destruction, therefore, clumps will not be able to enter the 
jet if Cc < Ww within the binary, which occurs at / > 0.1 
for the parameters given above. For relatively weak jets, 
say Lj < 10'^^ erg/s, these clumps can otherwise cross the 
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whole jet. For Lj ^ 10 erg/s, the clumps considered may 
enter the jet, but it would be difficult for them to escape 
(see Sect. 5.3). 



From the above discussion, we conclude that, under the 
jet-to- wind momentum flux ratios considered, if the wind is 
moderately inhomogeneous, say / < 0.1, clumps can effec- 
tively penetrate and mass-load the jet, with a penetration 
effectiveness depending on Re- As we show in what follows, 
clump penetration will have serious consequences on the 
jet stability and long-term collimation. We note that we 
have assumed that the properties of the clumps, and the 
mean separation between them, are constant all over the 
grid region in which they are located for simplicity. This 
approximation starts to fail for z > d, although most of the 
dynamical impact occurs below or around this distance. It 
is worthy to mention that for the purpose of this work, the 
assumed hydrodynamical nature of the jet could be loose 
to some extent, as long as the magnetic to kinetic pres- 
sure ratio in the jet is smaller, say, than ^ 0.1. Otherwise, 
the lateral magnetic pressure can effectively prevent clump 
penetration into the jet. Moreover, assuming that clumps 
managed to enter the jet, magnetic pressure could suppress 
the development of shocks, induce magnetic dissipation and 
non-thermal activity, and change strongly the cl ump dis- 
ru ption process, softening or potentiating it (e.g. Jones et 
ah] 1996 Shin et al. 2008). In any case, the flow evolution 
under a dynamically dominant magnetic field is out of the 
scope of this work. 



3. Simulations 

We have performed three numerical simulations using a 
finite-difference code named Ratpenat, which solves the 
equations of relativistic hydrodynamics in three dimen- 
sions, written in conservation form, using high-resolution- 
shock-capturing methods. Ratpenat was parallelized with a 
hybrid scheme with both parallel processes (MPI) and par 
allel threads (OpenMP) inside each process 



(see 



Perucho 

et al. 2010b). The simulations were performed in Mare 
Nostrum, at the Barcelona Supercomputing Centre (BSC) 
with 200 processors, each of them with a duration of 
1080 hours, amounting a total of 6.48 x 10^ computational 
hours. 



3.1. Simulation set-up 

The simulations are set-up with an overpressured jet sur- 
rounded by the stellar wind of the massive companion. We 
are thus implicitly assuming that the bow shock gener- 
ated by the jet when crossing this medium (PBKIO) is 
far enough and the cocoon dilute enough that the wind 
has occupied the space surrounding the jet itself]^ The 
physical size of the grid is, in (base) jet radius units, 
160 i?j X 160 i?j X 200 the last 200 i?j in the direction of 
propagation of the jet (z coordinate). The grid size in the x 
and y coordinates is divided in two regions, the inner 80 i?j 
around the jet axis having homogeneous resolution, and 
the outer 40 i?j in each direction being formed by cells with 
increasing size. The resolution in the homogeneous grid is 



^ The system we simulate is dynamic so the initial conditions 
are artificial. However, the long-term qualitative result is not 
affected by this, as the steady ]et is not in equilibrium but evolves 
governed by the interaction with the wind. 



4 cells/i?j at injectioifl with a total of 320 x 320 x 800 
cells. The extended grid is composed by 80 cells on each 
side, resulting in a box with 480 x 480 x 800 cells. In the 
simulations, _Rj = 2 x lO^'^ cm, so the physical size of the 
grid is (3.2 x 3.2 x 4) x lO^^ ^m. 

In previous works (PB08, PBKIO), the wind was simu- 
lated as homogeneous, but it was already suggested that a 
necessary improvement to those simulations should be the 
inclusion of inhomogeneities. We have done this by ran- 
domly adding Gaussian-shaped clumps in the side of the 
grid from where the stellar wind is injected. The initial 
number of clumps in the half of the homogeneous grid fac- 
ing the star, 1500, has been calculated to ensure that the 
mean wind mass-loss rate is « lO^^Afo/yr. The clump 
peak density is 10 times larger than the mean density 
(9 X 10~^^g/cm^), i.e., 9 x 10~^^g/cm^, and the minimum 
density between clumps has been fixed to 3 x 10~^^g/cm^. 
The whole wind region is set in pressure equilibrium (P^ = 
1.5 X 10"'^ erg/cm'^) with the higher density regions, and 
a velocity — 2 x 10^ cm/s radial from the star (PB08, 
PBKIO). Figure [2] shows cuts of the initial conditions in 
pressure and density along the axis of one of the simulated 
jets, and Fig.[3]shows transversal cuts at z ~ 2.2 x 10^^ cm. 



Log. rest- mass dais ily 





Fig. 3. Transversal cuts of the jet at half grid of axial ve- 
locity (cm/s) and rest-mass density (g/cm'^) for one of the 
initial models. 



We have simulated two jets with different powers, jet A 
with Lj = 3 X lO^*^ erg/s and jet B with Lj = 10^'^ erg/s (to 
be compared with jet 2 in PBKIO). The jets are set-up with 
a given opening angle but not in pressure equilibrium with 
the ambient, as this has homogeneous pressure whereas the 
jet pressure decreases with z. The injection point in the grid 
would be located at 2 x 10^^ cm from the compact object, 
this is, well within the binary system. The jets have both in- 
jection densities of 4.2 x 10~-^^g/cm^ and 1.4 x 10~-^''g/cm^, 
and a velocity uj = 10^° cm/s. Th ermal cooling terms, fol 



lowing the approximation used in jMyasnikov et al. (11998), 



have been added to th e code to account for the cooling in 
the clumps. Table [XT] summarizes the jet and wind param- 
eters used in the simulations. 

We note that to simulate continuous wind injection, the 
up-wind section of the ambient medium is replenished with 
clumps when a portion of its volume has been emptied of 
the original ones. This is done without any effect on the 
dynamics of the system, as the clumps are added in a re- 
gion far from the interaction between the jet and the first 
clumps. 



^ Note that the jet radius increases with 2, and so the effective 
resolution across the jet does. 
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Fig. 2. Axial cuts of the jet in the planes perpendicular (left) and parallel (right) to the jet-star plane of rest-mass density 
(top, g/cm'^), pressure (mid, erg/cm'^) and Mach number (bottom) for one of the initial models. 
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Table 1. Jet and wind parameters. 



A third simulation was performed to illustrate the evolu- 
tion of 3 clumps interacting with jet A (jet A'). The clumps 
were located at z = 3, 6 and 12 x 10^^ cm, with increasing 
distance to the jet in the x-direction such that the different 
clumps start their interaction with the jet subsequently. In 
Figure [4j the initial distribution of the clumps is shown in 
the axial density cut of the jet parallel to the star-jet plane. 



Logarithm of rest-mass density. Parallel to Jet-star plane. 




5.0X10" 1.0x10'^ 1.5X10'^ 2.0x10'^ 
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Fig. 4. Axial cut parallel to the star-jet plane of rest-mass 
density {g/cm^} for jet A' showing the initial location of 
the clumps. 



4. Results 

Jets initially expand due to their overpressure. The expan- 
sion is asymmetric, being faster in the direction of prop- 
agation of the wind (down- wind). The clumps that enter 
the jet close to its base are shocked by the expansion of 
the overpressured jet and eroded by its transversal velocity. 
The shock generated by this interaction propagates inside 
the jet, deforming its surface and triggering helical displace- 
ments of the jet flow, which add to the deviation of the jet in 
the direction of the wind further up. The material ablated 
from the clumps is dragged by the jet, circulates around its 
surface and generates a turbulent thick shear-layer where 
is mixed with material of the jet and accelerated in the 
direction of the latter. 

In the central region of the jet, a reconfinement shock 
is triggered by the irregular ram pressure of the clumpy 
wind. At this shock, the jet flow is decelerated, favoring 
clump penetration and clump-jet material mixing down- 
stream. At higher values of coordinate z, clumps can fully 
penetrate into the jet generating bow shocks that cover 
large fractions of the jet cross section and further decel- 
erate the flow. Thus, the nature of the jet is completely 
changed within the binary system. 
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4.1. Jet A 

In the case of jet A, the simulation reproduces ~ 1700 s of 
the jet-wind interaction. Due to the high overpressure of the 
jet at the injection point, the clumps are destroyed when 
they get close to the flow at the jet base, although much 
farther up they can penetrate inside the jet. Figures [5j [6] 
and [7] show two jet axial cuts, transversal and parallel to 
the star-jet plane, of rest-mass density, axial velocity and 
Mach number, respectively. The different images show the 
temporal evolution, with time increasing downwards. The 
initial phase of the interaction between the clumps and the 
jet is seen in the upper panels of each of the figures. The 
density maps show the generation of an asymmetric recon- 
finement shock (right panel) due to the ram-pressure of the 
clumpy wind, and the generation of a rarefied region in the 
boundaries of the jet produced by the initial expansion into 
the wind diluted regions. In this slower and more dilute rar- 
efaction region, clumps penetrate easily, deforming the jet 
surface (see Fig. [6]), and the up- wind jet expansion becomes 
dominated by the interaction with the clumps. 

At the adopted initial conditions, the jet cools down 
as it expands along the z axis, accelerating also slightly in 
the z-direction. This implies an increase of the Mach num- 
ber from the value at injection, M = 10, to M ~ 60. In 
the Mach number maps, we can see that, as the time goes 
by, the maximum Mach numbers in the jet downstream of 
the reconfinement shock get much smaller than the initial 
values. Deceleration due to mass entrainement, and heat- 
ing due to shocks, lead to the formation of transonic and 
subsonic regions inside the jet flow. 

The bottom panels of Figs. [5] [6] and [7] show the last 
snapshot of the simulation of jet A. The cuts transversal to 
the star-jet plane (left panels) reveal the jet expansion at 
small z, reconfinement and, once the jet has been deviated, 
the turbulent region where clump-jet mixing and jet decel- 
eration occur. The cuts parallel to the star-jet plane (right 
panels) show the effect that a clump entrained just after 
the reconfinement shock may have, triggering a deviation 
in the jet direction of more than 10°. As noted above, in 
the transversal cut to the star-jet plane of rest-mass density 
(Fig. [5] bottom left panel) , it can be seen that the clumps 
approaching the jet close to the injection point, where the 
overpressure is large, are destroyed before entering it. In 
this region, the figures also show complex structures of de- 
stroyed clumps mixed with jet material. 

Figure |8] shows cuts of rest-mass density, axial veloc- 
ity and Mach number transversal to the jet axis for jet 
A at the injection point, ^ 1/4, 1/2, and 3/4 of the grid 
length (z = 2 X 10", 1.2, 2.2, and 3.2 x lO^^ ^m). The first 
cut shows the disruption of the clumps as soon as they 
get close to the jet, whereas the jet Mach number and 
velocity keep their initial values. At one quarter of the 
grid {z ~ 1.2 X 10"'^^ cm), the jet presents an elongated 
structure in the down-wind direction. In the up-wind di- 
rection, rarefied jet material left after the initial jet ex- 
pansion forms filaments between the clumps. At this point, 
some regions of the jet section keep a large fraction of the 
original Mach number and axial velocity, but some oth- 
ers already show signs of the action of the reconfinement 
shock. At z ~ 2.2 x 10^^ cm, the jet flow has already gone 
through this shock, and the wind thrust has managed to 
deviate its original direction. Although the flow has reex- 
panded and partially recovered the strongly supersonic na- 



ture, the impact of the reconfinement shock has caused a 
noticeable decrease in velocity from the original values. At 
z ~ 3.2 X 10^^ cm the jet shows a large radial displace- 
ment and irregular morphology. Overall, the jet section is 
clearly displaced from the center of the grid. Although the 
jet keeps a high axial velocity at this position, the density 
cuts show an increase of more than one order of magnitude 
with respect to the original value (~ 3.2 x 10~^^ g/cm~'^) 
in some regions of the jet at this position. 

4.2. Jet B 

This simulation reproduces the jet-wind interaction during 
t ~ 2000 seconds for jet B, with Lj = lO^"^ erg/s. The 
same qualitative explanation for the evolution of jet A ap- 
plies also to jet B. However, the overpressure of the jet 
is larger, and the reconfinement shock, produced by the 
global clumpy wind impact, appears further downstream 
than in jet A. There are interesting differences with re- 
spect to jet A due to the entrainement of a clump in jet B 
at z ~ 2.2 X 10^^ cm, which is favored by the strong de- 
crease of the inertia of the jet downstream of the recon- 
finement shock. This clump generates a strong bow shock 
that eventually covers the whole jet section. Behind the 
bow shock the jet is decelerated and a turbulent layer is 
formed. Once the clump is destroyed, the jet may recover 
its previous configuration, but the constant injection of fur- 
ther inhomogeneities in the wind makes it unlikely. This 
event is very illustrative of the clumpy wind disruptive 
effects. Figures (Oj [lO] and [TT] show the evolution of jet B 
with time through panels of rest-mass density, axial veloc- 
ity and Mach number, respectively. The large overpressure 
of the jet allows some filaments of jet material to propa- 
gate between clumps in the up-wind direction. Again, the 
jet becomes subsonic downstream of the region where the 
reconfinement shock and the bow shocks generated by the 
clumps occur (z ~ 2.4 x 10^^ cm). 

Figure [12] shows transversal cuts of axial velocity, Mach 
number and rest-mass density at the last snapshot of the 
simulation at the jet injection position, z ~ 2 x 10^^ cm, and 
at z = 1.2, 2.2, and 3.2 x 10^^ cm. At the injection point, the 
image obtained is very similar to that from jet A (see upper 
panels in Fig. [s]). At z = 1.2 x 10^^ cm and 2.2 x 10^^ cm, 
the jet is still well collimatcd with large velocity and Mach 
number. In these cuts, filaments of material travelling in the 
up-wind direction can be observed. The last set of images 
(bottom row), at z = 3.2 x 10^^ cm, shows the effect of the 
entrained clumps. In the left panel, the traces of the clumps 
can be seen as low velocity regions, resulting also in tran- 
sonic or even subsonic velocities as indicated by the Mach 
number maps. The initial density of the jet at this point 
was p ~ 5 X 10~^^g/cm~^ at the start of the simulation, 
whereas the right panel shows densities one order of magni- 
tude larger within the jet, indicating efficient assimilation 
of wind material. 

4.3. Jet A 

The previous simulations show very complex structures and 
dynamics. In order to have a clearer idea of the processes 
taking place, we repeated the simulation of jet A with only 
three clumps located initially at z = 0.3, 0.6 and 1.2 x 
10^^ cm, as shown in Fig. |4j The simulation was stopped 
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Fig. 5. Axial cuts of rest-mass density (g/cm^) at different instants for jet A. Top panels stand for t — 300 s, second row 
for t — 1220 s, t = 1570 s the third one, and bottom one for the last snapshot of the simulation, at i = 1700 s. 



after t ~ 1000 seconds. Figures [13] and [M] show the time 
evolution of the interaction of the clumps with the jet. The 
first clump is destroyed in the jet surface, with little dy- 
namical impact on the jet. Otherwise, the second clump 
triggers a shock that propagates inside the jet, although 
the clump itself is also destroyed before entering into the 
jet. Between the first two clumps, jet material shocked by 
the second clump rises due to the local pressure gradients 
in the up-wind direction. The third clump shows the same 
behavior as observed in the simulation of jet B, being able 
to enter into the jet before being destroyed, and generating 
a strong bow shock that causes jet expansion, deceleration 
and significant mass-load. 

Figure [15] shows transversal cuts of axial velocity, Mach 
number and density (as in Figs. Js] and 12). In this case, the 
cuts are done at z ~ 0.25, 0.75, T, 1.5 and 1.75 x 10"^^ cm. 
The structures are much clearer in these plots due to 
the small number of clumps. At z = 0.75 x 10^^ cm we 
see the disrupted rests of the first clump (originally at 
z — 0.3 X 10^^ cm, see bottom panel of Fig. 13). Initially, 



this clump crossed the shock driven by the jet in the di- 
luted medium, and was disrupted at the jet surface. This 
interaction already implies some loss of initial jet thrust, 
which is transferred to the clump remains that become part 



of the jet flow. The shock that this interaction generates is 
observed as a bow-like structure inside the jet inner bound- 
ary in the Mach number and density panels at higher values 
of coordinate z. The next set of cuts, at 10^^ cm, shows a 
snapshot of the disruption of the second clump. This pro- 
cess is here at an earlier stage than that of the first clump. 
At z = 1.5 and 1.75 x 10^^ cm, one can see the clear effect of 
the clump located at 1.4 x 10^^ cm, which generates a bow- 
like structure entering the jet and a cometary tail, crossed 
by the transversal cuts, of jet dragged material (see bottom 
panels in Figs. 13 and 14 ). The other bow-like waves, visible 



mainly at x < in the images, are the bow shocks gener- 
ated by the interaction with the clumps upstream of this 
position. At z = 1.75 x 10^^ cm, the jet is dominated by the 
bow-like structures triggered by the interactions upstream, 
and the cometary-like tail from the third clump is seen at 
X > region. The faster regions of the jet are those that 
have not been crossed by any of these shocks, i.e., those that 
lie at the smaller values of x. The jet keeps a large velocity 
in those particular regions, but it is in general efficiently 
mass-loaded and decelerated. 

Overall, the picture is similar to that shown in simu- 
lations of jets A and B, showing that only a few clumps 
generated by inhomogeneities in the injection of the wind 
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Fig. 6. Axial cuts of axial velocity (cm/s) at different instants for jet A. Top panels stand for t — 300 s, second row for 
t — 1220 s, third one for t = 1570 s, and bottom one for the last snapshot of the simulation, at t — 1700 s. 



could produce already significant changes in the properties 
of the jet flow. 



5. Discussion 

5.1. Jet propagation 

The simulated jets are strongly deviated from their original 
direction and show disrupted morphologies at the end of the 
calculations. Initially, the first interactions at small values 
of z trigger helical patterns that, in the absence of further 
perturbations, could couple to Kelvin-Helmholtz unstable 
modes. However, the reconfinement shock triggered by the 
ram pressure of the clumpy wind on the jet plus individ- 
ual interactions with multiple clouds downstream of this 
shock, generate non-linear structures that are most impor- 
tant for the jet evolution. The position of the reconfinement 
shock is located within the binary region, as predicted by 
PB08. There, the jet is decelerated and heated, thus be- 
coming even weaker relative to the wind thrust, as also 
shown in PB08 and PBKIO. In the present case, the clumps 
from the up-wind region can penetrate in the jet already at 
z ~ 1 — 2 X 10^^ cm, generating additional bow-like shocks 
and decelerating even more the jet flow. The deviation of 
the jet favors the presence of a region, on the down-wind 



side of the jet, to which, by the end of simulations, the 
clump bow shocks and jet mixing have still not reached. In 
this region, the flow keeps a relatively large axial velocity 
in this direction. However, both shocks and clump-mixed 
jet material should eventually fill the whole jet at larger z, 



so the description just presented (see Figs [9 10 and 11) 
can be extended to the whole jet cross section. As a result, 
farther downstream the jet becomes mass-loaded, slow and 
transonic. Although a similar result was already obtained in 
the homogeneous wind case (PB08, PBKIO), wind clump- 
ing significantly potentiates jet destruction. 

When comparing jet B in this work with jet 2 in PBKIO, 
it is remarkable that the presence of a forward shock and 
the absence of clumps in the latter have a strong influ- 
ence on the long-term evolution of the jet. The jet-driven 
forward shock in the wind and the corresponding cocoon 
(see PBKIO) keep the wind at some distance from the jet. 
If a homogeneous wind were already in contact with jet 2 
in PBKIO, a smooth shear layer would form. Otherwise, 
the presence of clumps in jet B enhances the thrust lo- 
cally in the wind, increasing mass, momentum and en- 
ergy exchange. Jet 2 in PBKIO thus shows a larger de- 
gree of collimation and is just slightly deviated from its 
original direction of propagation, contrary to what we ob- 
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Fig. 7. Axial cuts of Mach number at different instants for jet A. Top panels stand for t = 300 s, second row for t = 1220 s, 
third one for t = 1570 s, and bottom one for the last snapshot of the simulation, at ^ = 1700 s. The color scale for this 
parameter has been saturated at M = 5, thereby enhancing the contrast inside the jet and thus the location of transonic 
flows resulting from shocks. 



serve in the case of jet B in this paper. From this work 
we can conclude that jet luminosities Lj > 10'^'' M-gerg/s 
(M_6 = M/lO-s Mq/yt) are needed if the jet is not to be 
destroyed when crossing the binary system. We notice that 
these dynamical arguments favor Lj > 10^^ and 10'^® erg/s 
in the HMMQ Cygnus X-1 and Cygnus X-3, respectively. 
The cocoon/clumpy wind case deserves few words. This sit- 
uation takes place when the forward shock is well within the 
binary system, and the cocoon pressure is high, preventing 
clumps from entering the cocoon. When the forward shock 
has reached the outskirts of the binary, the cocoon pressure 
drops quickly (PBKIO), allowing wind clumps to penetrate 
into the cocoon and reach the jet, and eventually dissipate 
the cocoon away. 

5.2. Clump evolution 

In jets A and B, those clumps reaching the jet relatively 
close to its base are destroyed just by jet expansion or ero- 
sion. These interactions trigger a shock wave that propa- 
gates inside the jet, and when interactions are frequent such 
a wave forces the strongly disruptive asymmetric recoUima- 



tion shock. This phenomenon is illustrated in Fig. 13), in 



which the first clump is completely destroyed by jet ex- 
pansion, whereas the second one, even when it does not 
fully penetrate into the jet, triggers a shock strong that 
propagates all through the latter lasting for the whole sim- 
ulation (several t^; see Sect. [2]). The clump at the high- 
est z in jet A' (1.4 x 10^^ cm) is shocked but still not 
significantly disrupted by instabilities after few td- Later, 
this clump could be destroyed or may eventually escape 
the jet, although jet bending in the down-wind direction 
makes the latter unlikely. When clumps are disrupted in- 
side the jet, all the clump mass is entrained by the flow. 
The level of mass loading can be easily estimated from the 
amount of clumps entering into the jet per time unit: A^cj ^ 
(77/47r)(3M/47ri?3pJ « 0.02 clump/s or « 5x10^^ g/s. This 
is 3 times more mass flux than in the jet. Therefore, jet 
deceleration due to mass loading can be very efficient, as 
is most clearly seen in the simulation results for jet B. The 
implications also apply to faster and lighter, but equally 
powerful jets. Even if a strongly relativistic flow is injected 
at the jet base, for similar jet (relativistic) ram pressures 
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Fig. 8. Transversal cuts of rest-mass density (g/cm^), axial velocity (cm/s) and Mach number for the last snapshot 
{t = 1700 s) of jet A, at z = 2 X 10^^, 1.2, 2.2, and 3.2 x 10^^ cm. The color scale of the Mach number has been saturated 
at M 5. 

clump penetration into, and mixing with, the jet will also 5.3. Radiation considerations 
occur. This will mass-load and brake the jet very quickly. 



The qualitative dynamical scenario presented in ABR09 
is validated by our numerical work in a semi-quantitative 
way. Remarkably, as hinted in ABR09, the clump interac- 
tion with many clumps will reduce the jet kinetic luminosity 
and also its ram pressure. In this way, clumps can actually 
keep their integrity longer, and penetrate farther inside the 
jet, even for moderate /-values. 



In ABR09, the radiation produced by a clump inside a 
HMMQ jet was calculated. Given the compactness of the 
considered region, i.e. the bow shock formed around the 
clump, radio emission was negligible. The luminosities of 
synchrotron X-ray, and (stellar photon) IC gamma-rays, 
dominated the non-thermal output, with their values anti- 
correlating depending on the magnetic field. Based on that 
work, and extending the study to the multi-clump /jet in- 



9 



Perucho & Bosch-Ramon: 3D simulations of jets in clumpy winds 



Logarithm of rest-mass density. Transversal to jet-slar plane. 




2x10 3x10 
z (cm) 



Logarithm of rest-mass density. Transversal to jet-slar plane, 




2x10 3x10 
z (cm) 



Logarithm o! rest-mass density. Transversal to jet-star plane. 




2x10 3x10 
z (cm) 



Logarithm of rest-mass density. Transversal to jet-star plane. 




2x10 3x10 
z (cm) 



Logarithm of rest-mass density. Parallel to jet-star plane. 

"1 




Logarithm of rest-mass density. Parallel to jet-slar plane. 



4.73e-13 












1 7.19e-ie 




HI 




1 




1.09e-18 




-5>t10" 1 



1x10 2x10" 3x10 4x10 

z (cm) 



Logarithm of rest-mass density. Parallel to jet-star plane. 




Logarithm of rest-mass density. Parallel to jet-star plane. 




Fig. 9. Axial cuts of rest-mass density {g/cm^^) at different instants for jet B. Top panels stand for t = 470 s, second row 
for t — 1640 s, third one for t — 1770 s and bottom one for the last snapshot of the simulation, a.t t = 1970 s. 



teraction case simulated here, we qualitatively discuss in 
what follows the expected emission. 

The typical number of clumps inside the jet, at the bi- 
nary scales, can be estimated as Ncj ^ few times td A^cj, i-e. 
> 10 for the wind and jet properties adopted in this work 
(recall that the real clump lifetime is longer than t^). Since 
td ~ 100s < i?j/ww ^ 1000 s, clump destruction will likely 
limit iVj instead of jet-crossing escape. Although the jet 
can be quite disrupted and mass-loaded at the borders of 
the system, one can derive an estimate for the non-thermal 
emission, which depending on B'^/Sttu^, (m* = L^/ATr(Pc) 
will be released either through synchrotron (keV-MeV) 
or IC (GeV-TeV). We note that, for values well below 
equipartition with the non-thermal electron energy density, 
synchrotron self-Compton may also contribute to the IC 
component, although we neglect this component here for 
analysis simplicity. We also do not consider non-radiative 
losses or hadronic processes, but refer to ABR09. 

Assuming a 10% efficiency for the kinetic-to-non- 
thermal energy transfer in clump bow shocks, and account- 
ing for the clump-jet covering fraction ^cj {Rc/vd)^7 
the high-energy emission can reach luminosities of ~ 
0.1 Ncj ^cj Lj ^ 0.01 Lj. This value may lead to a source 
that may be detectable in GeV {B"^ /Attu^, < 1) and even 



in TeV {B'^/4ttu^ ^ 1), sho wing persistent plu s fluctu- 



ating variability components (Owocki et al. 2009). Strong 



TeV photon absorption may be avoided since most of the 
emission would come from the bor ders of the binary (e.g., 
Bosch- Ramon fc Khangulyan|2009 ). For B^ /'hiu^ > 1, syn- 
chrotron soft gamma-rays would dominate the non-thermal 
output. Thermal radiation, peaking at several keV, from 
the shocked clumps may be also present ( ABR09) . but it is 
minor for the adopted clump properties. 

It is noteworthy that for a steep clump-size distribution 
(see Sect. [2]), i.e. with most of the clump-jet interaction 
occurring in a quasi-homogeneous wind regime, still several 
big clumps could enter the jet. This would likely lead to a 
lower level of persistent emission, but time to time chance 
would bring few big and dense clumps together inside the 
jet, rendering sudden and bright events with durations of 
a few 1000(i?c/10"cm) (cc/lO^cm)-! s (ABR09). These 
statistically expected violent events could be related to a 
temporary complete destruction of the jet for large enough 
clumps. 

Although the jets of microquasars present a com plex 



phenomenology of their own ( Fend er et al.[2004 2009 1 , on 
top of that it is expected additional phenomena linked to 
the jet-wind interaction in HMMQ, as this and previous 
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Fig. 10. Axial cuts of axial velocity (cm/s) at different instants for jet B. Top panels stand for t — 470 s, second row for 
t — 1640 s, third one for t = 1770 s and bottom one for the last snapshot of the simulation, a,t t = 1970 s. 



works (PB08, PBKIO) show. Since the jet power is strongly 
linked to the mass-loss rate, HMMQ jets detectable in X- 
rays or gamma-rays will be probably found in systems with 
moderate-to-strong winds. Therefore, unless particle accel- 
eration is negligible in clump bow shocks, HMMQ phe- 
nomenology at high-energies must be strongly affected by 
them, showing wind-related strong variability at high ener- 
gies unless the wind is quite homogeneous or clumps with 
i?c > 10^'^ cm are completely missing. Concerning (pow- 
erful) transient ejections, usually associated to X-ray state 
transitions, we note that such an ejections will require some 
time to form. If this takes hours, the wind will have time to 
surround the transient jet. Then, the clump impact will be 
as described here unless the jet is too powerful. If power- 
ful blobs would appear as discrete even at the scales of the 
binary, they may have too much inertia to be significantly 
affected by the wind. 

The dynamical impact of the stellar wind on the jet, 
enhanced by dumpiness, should not be neglected when in- 
terpreting radio emission from HMMQ. Even if jets escape 
from disruption, the enhanced jet entropy and bending, 
even by small angles 10°), could have observational con- 
sequences. The reason is that bending pushes jet material 
farther from the orbital plane axis, increasing the strength 



of the Coriolis force exerted by the stellar wind and the 
orbital motion. The farther away from this axis the jet 
reaches, the stronger this force gets, enhancing jet heating, 
turbulence and bending. Although quantitative predictions 
call for a detailed study, this region may be observation- 
ally probed using VLBI techniques. Far enough from the 
binary, once the jet has become too wide to be affected by 
the orbital motion, a coUimated supersonic flow may form 
again. If the energy and momentum fluxes kept enough 
anisotropy after crossing the system and suffer from or- 
bital motion, given a negligible external pressure a jet-like 
structure could form again and propag ate unstopped up to 
pc scales, terminat ing in th e ISM (e.g. Bordas et al.||2009 



'B osch-Ramon et ar.||2011t |Yoon et ahpOlip . Despite the 
different outflow geometry, similar phenomena are also ex- 
pect ed in high-mass binaries hostin g non-accreting pulsars 
(see Bosch-Ramon fc Barkov|[2011 ). 
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Fig. 13. Axial cuts of rest-mass density (g/cm'^) at different instants for jet A'. Top panels stand for t = 170 s, second 
row for t = 320 s, third one for t = 780 s, and bottom one for the last snapshot of the simulation, a,t t = 970 s. 
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Fig. 14. Axial cuts of axial velocity (cm/s) at different instants for jet A'. Top panels stand for t — 170 s, second row 
for t — 320 s, third one for t = 780 s, and bottom one for the last snapshot of the simulation, a.t t = 970 s. 
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Fig. 15. Transversal cuts of rest-mass density (g/cm^), axial velocity (cm/s) and Mach number for the last snapshot 
{t = 970 s) of jet A' at z ~ 0.25, 0.75, 1, 1.5 and 1.75 x 10^^ cm. The color scale of the Mach number has been saturated 
at M = 5. 
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